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Abstract 

Employee stock options (ESOs) are American-style call options that can be terminated early 
due to employment shock. This paper studies an ESO valuation framework that accounts for job 
termination risk and jumps in the company stock price. Under general Levy stock price dynamics, 
we show that a higher job termination risk induces the ESO holder to voluntarily accelerate 
exercise, which in turn reduces the cost to the company. The holder’s optimal exercise boundary 
and ESO cost are determined by solving an inhomogeneous partial integro-differential variational 
inequality (PIDVI). We apply Fourier transform to simplify the variational inequality and develop 
accurate numerical methods. Furthermore, when the stock price follows a geometric Brownian 
motion, we provide closed-form formulas for both the vested and unvested perpetual ESOs. Our 
model is also applied to evaluate the probabilities of understating ESO expenses and contract 
termination. 
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1 Introduction 


Employee stock options (ESOs) are an integral part of executive compensation in the United States. 
The prim ary objecti v e is t o align the interests between the executives and the firm. According to 
Frvdman and .Tenter ( 201C)h . stock options account for 25% of the total compensation package of 
CEOs in 2008. In Table 1, we see that the percentage of S&P 500 companies granting ESOs remains 
above 73% in 2011, despite a decline from 93% in 2000. The cost of these options can potentially 
be very burdensome to the shareholders. In view of this, the Financial Accounting Standards Board 
(FASB) has passed regulations to require firms to estimate and report the granting cost of ESOs. 

Typically, ESOs are early-exercisable long-dated call options written on the company stock. To 
maintain the incentive effect, the firm usually imposes a vesting period that prohibits the employee 
from exercising the option. During the vesting period, the employee’s departure from the firm will 
result in forfeiture of the option (i.e. it becomes worthless). After the vesting period, when the 
employee leaves the firm, the ESO will expire though the employee can choose to exercise if the 
option is in the monej0. Table [1] summarizes the average vesting period and average maturity of 
ESOs granted by S&P 500 companies over 2000-2011. As we can see, the vesting period has been 
consistently close to 2 years, while the average maturity has decreased from 9-10 years to 8 years over 
a decade. 

The key to ESO valuation involves modeling the employee’s voluntary exercise strategy as well 
as job termination time, especially since the option is typically terminated prior to the contractual 
expiration date. Moreover, the possibility of future employment shock can influence the employee’s 
decision to exercise now or later. In fact, the FASB guideline]! also recommends that any reason¬ 
able ESO valuation model incorporate “the effects of employees’ expected exercise and post-vesting 
employment termination behavior.” 

In this paper, we study a valuation framework that incorporates the common ESO features of 
vesting period, early exercise and job termination risk, while allowing for different price dynamics 
with jumps. Specifically, we model the arrival of the employment shock by an exogenous jump 
process, and formulate the American-style ESO as an optimal stopping problem with possible forced 
exercise prior to expiration date. Our valuation problem is studied under a class of exponential 
Levy pric e processes, rath e r than limiting to the GBM fr amework commonly fou nd in the literature 


(see e.g. Hull and White ( 20041 ) : Cvitanic et al.1 ( 20081 ) : Carpenter et al. ( 201fll )). Under different 


job termination intensity assumptions (constant or stochastic), we analyze the corresponding free 
boundary problems in terms of an inhomogeneous partial integro-differential variational inequality 
(PIDVI), and discuss the computational methods to solve for the option value as well as the optimal 
exercise strategy. Analytically and numerically, we find that with higher job termination risk it is 
optimal for the holder to voluntarily accelerate ESO exercise. For risk analysis, we also apply our 
numerical schemes to calculate the probability of cost exceedance and the probability of contract 
termination under various scenarios. 

In existing literature, there are several major approaches to model the exercise timing of ESOs. 
The first approach assigns an ad hoc exercise boundary for the ESO holder. In other words, the ESO 
exercise is characterized by the first passage time of underlying stock price to a pre-specified level, 
which is independent of the employee’s job termination rate or other model parameters. The main 
advantage of this approach is the availability of closed-form formu las for the ES O valu e. A ddition ally, 
one can also incorporate a random job termination time, as in Hull and White ( 2004 1: Cvitanic et al. 


1 See, for example, “A Detailed Overview of Employee Ownership Plan Alternatives” by The National Center for 
Employee Ownership (http://www.nceo.org). 

2 See Sect. A.16, FASB Statement 123R (revised 2004), available on http://www.fasb.org/summary/stsuml23r.shtml. 
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Year 

% of S&P500 companies 
granting ESOs 

Avg. vesting period 
(years) 

Avg. Maturity 
(years) 

2000 

92.98% 

2.00 

9.24 

2001 

94.35% 

2.22 

9.28 

2002 

93.56% 

2.18 

9.53 

2003 

89.59% 

2.18 

10.17 

2004 

88.09% 

2.03 

8.66 

2005 

75.34% 

2.16 

8.61 

2006 

81.30% 

2.12 

7.86 

2007 

77.33% 

2.18 

8.14 

2008 

76.31% 

2.38 

7.35 

2009 

75.30% 

2.16 

8.41 

2010 

72.69% 

2.32 

8.71 

2011 

73.39% 

2.16 

8.07 


Table 1: Summary of ESO compensation during 2000-2011. Source: Thomson Reuters Insider Trading database 
and Compustat ExecuComp database. 


( 20081 ). but it would not directly affect the voluntary exercise boundary. 

To incorporate the employee’s risk attitude, the utility maximization approach derives the em¬ 
ployee’s exercise strategy from a p ortfolio optimization problem that accounts for the employee’s risk 
avers i on and hedging re s trictions (IHuddartl. 1 1994; IChance and Yana . 120051 : iGrasselli and Hendersonl . 


20091 : iLeung and Sircad. [2009|; I Carp enter et all l201flh . The utility-maximizing exercise boundary 
represents the optimal voluntary exercise timing for the risk-averse employee. From the firm’s per¬ 
spective, the ESO cost is computed by risk-neutral expectation assuming the option will be exercised 
at the utility-maximizing boundary or job termination time, whichever is earlier. 

On the other hand, the intensity-based approach does not distinguish voluntary and involuntary 
exercise, and model the contract termination by a single r andom time, characterize d by t he first 
arrival time of an exogenous jump process. In the literature, Jennergren and Na sl und ( 19931 ) model 
the exercise time using an exogenous Poisson process, and Carr and Linetskv (j2QQoh propose an 
intensity-based model for ESO valuation where both job termination a nd volunta r y exe rcise intensities 
are functions of the company stock price. In the theoretical study by ISzimaverl ( 20041 ). the employee 
is allowed to optimally select a voluntary exercise time, while the sudden departure is represented 
by an exogenous Cox process. In practice, the job termination intensity specification depends on the 
firm’s history and estimation methodology. For empirical stud i es on the ea r ly exe rcise patterns of 
ESOs, we refer to Huddart and Lang ( 19961 ) : Marquardt ( 20021 ): Bettis et al. (2005). 

We consider the ESO as an American call option with a possible sudden forfeiture or forced 
exercise due to job termination. In contrast to many existing ESO models, we study a versat ile 
valuat ion f ramework t hat is compatible for a wide class of Levy price pr ocesses, including iMerton 
( 19761 ) and Kou ( 2002) jum p diffusions, as well as Variance Gamma (VG) ( Madan et al. ( 19981 )) and 
CGMY ( Carr et al. ( 20021 )) models, in addition to GBM that is commonly found in the literature 
including those cited above. This allows us to study the combined effect of jumps in stock price and 
job termination intensity on the ESO value and optimal exercise boundary. Moreover, our model 
provides the end-user, i.e. the firm, the flexibility in choosing the appropriate price model (within 
a general Levy class) for the underlying stock. Currently, the FASB permits the use of the Black- 
Scholes formula with the ESO’s contractual term replaced by its average life, as well as a number of 
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variation^. In this regard, our paper offers an alternative valuation approach that accounts for both 
voluntary exercise and job termination risk along with various choices of models for the company 
stock price. 

The ESO valuation proble m can be con si dered as an exten s ion to the pricing of Ameri c an op ¬ 
tions under Levy proce sses; see Pham ( 1997 1. Hirsa and Madanl ( 20041 ) . iBavraktar and Xind ( 20091 ) . 
Lord et al. ( 20081 ) . and Lamberton and Mikou 1 20081 ) . among others. For vested ESOs, the job ter¬ 
mination arrival forces the employee to exercise immediately. This leads to an optimal stopping 
problem with forced exit. In terms of the variational inequality for the option value, this restriction 
gives rise to an inhomogeneous term that depends on both the option payoff and job termination 
intensity. When the stock price follows a geometric Brownian motion, we provide the closed-form 
formulas for both the vested and unvested perpetual ESO costs. The optimal exercise threshold can 
be determined uniquely from a polynomial equation, and it admits an explicit expression in the case 
without job termination risk. 

In order to compute the ESO value and exercise boundary, we apply the Fourier Stepping Timing 
(FST) approach, whereby the associated inhomogeneous PID VI is sim plifie d by F ourier transform and 
the optimal exercise price is determined in each time step. Jackson et al. ( 20081 ) apply this approach 
to price European, American, and barrier options under a number of Levy models. For our ESO 
valuation problem, the structure of the inhomogeneous PIDVI varies under different job termination 
intensity specifications. In particular, if the intensity is affine in the (log) stock price, then the 
inhomogeneous PIDE for the ESO in the continuation region can be simplified to an inhomogeneous 
PDE. In the constant intensity case, we further reduce the associated PIDE into an ODE. These 
observations lead to several efficient numerical algorithms for valuation. In all these cases, we compare 
with the numerical results from an implicit-explicit finite difference method for valuing the ESO. To 
this end, we refer to the fini te difference methods for pr ic ing American op t ions under Levy or jump 
diffusi on m odels, includin g Cont and Voltchkova ( 2003l l: d’Halluin et al. ( 20031 ) ; Hirsa and Madan 
( 20041 ) . and Forsyth et al. ( 20071 ) . 

Firms typically expense the granted ESOs according to a fixed schedule, such as quarterly or 
annually, but the cost of an ESO changes over time depending on the stock price movement. From 
the firm’s perspective, a rise in the ESO value implies a higher expected cost of compensation as 
compared to the initially reported value. On the other hand, existing ESOs can be either exercised 
voluntarily by the employee, or terminated due to employment termination. This motivates us to 
study (i) the probability that the ESO cost will exceed a given level in the future, and (ii) the contract 
termination probability. The ESO cost exceedance probability bears similarity to the loss probability 
used in classical value-at-risk calculation. The contract termination probability sheds light on the 
likelihood that the firm will have to pay the employee over a future horizon, from a week to a few 
years. We apply Fourier transform based methods to compute these probabilities under constant 
job ter mi nati on i ntensit y, and we show the connection between our approach and that developed by 
Carr and Madan ( 19991 ) in the computation of these probabilities. 

The rest of the paper is structured as follows. In Section [2l we formulate the ESO valuation 
framework under Levy price dynamics. In Section[3l we discuss the Fourier transform based numerical 
methods for ESO valuation. In Section we discuss the valuation problem under stochastic job 
termination intensity. In Section [5j we analyze and evaluate the probability of ESO cost exceedance 
and the probability of contract termination. In Section [6l we provide closed-form formulas for the 
vested and unvested perpetual ESO costs under the lognormal stock price model. Section [7] concludes 
the paper. 


3 See Sect. A.25, FASB Statement 123R (revised 2004). 
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2 Model Formulation 


In the background, we fix a probability space (fi,(/, P) satisfying the usual conditions of right conti¬ 
nuity and completeness, where P is the historical probability m easu re. Let {X t )t >o be a Levy process, 
which admits the well-known Levy-Ito decomposition (S atol . 119991 . p.119): 


X t = yt + crB t + X[ + lim Y t e , 

e\0 


where B is a standard Brownian motion under 


X 0 — 0, 

and the jump terms are given by 


X[ = 


x e t = 


f 

J\v 


>l,sG[0,t] 


yj(dy,ds ), 


/ y(J(dy,ds) - v(dy)ds) = yJ(dy,ds). 

J e<|y|<l,se[0,t] e< |?/| <l,se [0,t] 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


The characteristic triplet (/r, cr 2 ,^) of X consists of the constant drift y and volatility a, along with 
the Levy measure v. In (12.2p and (12.31) . the Poisson random measure J(dy,ds ) counts the number of 
jumps of size y occurring at time s, and J is the associated compensator. 

The characteristic exponent T (u) of X is given by the Levy-Khintchine formula (jSatol . 1999, 
p.119): 

T(w) = iyu - \(t 2 u 2 + [ (e l “ y - 1 - iyul { \ y \ <1} ) v(dy), u <E C. (2.4) 

z J K 

With this, the characteristic function of X t is 4>x t (u) = We denote F A = (T^)t >0 as the 

filtration generated by X. 

In Table El we summarize the Levy densities and characteristic exponents f or sever al well -known 
Levy models that will be used in this paper. In the GBM model ([Black and Scholes . 1973l h the 
Levy density is absent as the stock price has no jumps. The Merton ( 19761 ) jump diffusion model 
features normally distributed jump sizes, while the[Kouj (120021 ) model assumes a doubl e expo nential 


distribution for the jump sizes. In con trast, under the Variance Gamma (VG) model (i Mad an et al 


( 1998l ll and CGMY model ( Carr et al. ( 2002 1). the stock price follows a pure jump process with 
infinite activity. 


Model 

Levy Density v(dy) 

Characteristic Exponent T(u;) 

GBM 

N/A 

iyu — 

Merton 

a c -b({v-d)/v) 2 

V2tto 2 

iyu <j2 “" + a(e lfluJ " 2 1) 

Kou 

VG 

a(prj + e- yri +l{ y>0 }+ (1 - J»)»7- e“l| w<0} ) 
l e biy-b 2 \y\ 

KF+w( e l ' 7ll {y< 0}+e / l{ 2/ >0}) 

i yu 2 +a(i_ iu)/?7+ + X ) 

l log(l i~yKu + ° 2 ™ 2 ) 

CGMY 

CT(- Y) [{M - iu) Y - M y + (G + iu) Y — G y ] 


Table 2: Levy densities and characteristic exponents for some well-known Levy processes. In all these models, 
y and a are the drift and volatility of the Brownian Motion under P. For the Merton and Kou models, a 
denotes the jump intensity. In the Merton model, y and <r 2 are the mean and variance of the IID normally 
distributed jumps. In the Kou model, p (resp. (1— p)) is the probability of positive (resp. negative) exponential 
jumps. In the VG model, b\ = y/d 2 , 62 = \Jy 2 + 2<f 2 /n/d 2 , and in the CGMY model, C,G,M > 0, Y <2. 

The company stock price is modeled by an exponential Levy process 

S t = S 0 e Xt , t > 0, 
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with constant initial stock price So > 0. In addition, we assume positive constant interest rate r and 
non-negative dividend rate q. For ESO valuation, we work with a risk-neutral pricing measure (Q> 
such that 

E^{ e Xl } = e r ~ q <0> 'F(— i) = r — q, 

where T is given in (12.41) with n replaced by 

<7^ f 

fi = r - q - — - / (e y -1- yl{\ y \<i}) i>(dy), 

A Jr 

and D is the Levy measure under Q. 

2.1 Payoff Structure 

Figure [T| illustrates the payoff structure of an ESO with strike K, vesting period of t v years, and 
expiration date T. During the vesting period, the ESO is not exercisable and is forfeited if the 
employee leaves the firm. As soon as the option is vested (at or after time t v ), the employee can exercise 
the option at any time prior to the expiration date T, but will also be forced to exercise immediately 
upon job termination. We model the employee’s job termination time r A by an exponential random 
variable with rate parameter A > 0, and assume that r A and X are independent. The case of 
stochastic job termination will be discussed in Section 0J 


(2.5) 

( 2 . 6 ) 



Figure 1: ESO payoff structure. From bottom to top: (i) The employee leaves the firm during the vesting 
period, resulting in forfeiture of the ESO. (ii) The employee is forced to exercise the vested ESO early due to 
job termination, (iii) The stock price path jumps across the exercise boundary after vesting, so the employee 
exercises the ESO immediately, (iv) The stock price reaches the exercise boundary (without jump) after 
vesting, and the option is exercised there, (v) The employee exercises the ESO at the end of vesting. Along 
each stock price path, the vertical line segments depict jumps in the stock price. 
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(2.7) 


2.2 ESO Cost 

The value of a vested ESO at time t € [ t v ,T] is given 1 

C(t,x) = sup Eg / e -r(rAr A -t)( jS r oe Jf T ArA -K)+ 
r&Tt,r ’ l 

= sup Eg / e - (r+A)(r - t) (S 0 e^ - K) + + 

T&T t ,T ’ l 

where 7 t,T is the set of F A " stopping times taking values in [t, T], and Eg,,{-} denotes the conditional 
expectation with Xt = x. In other words, after the vesting period, the employee faces an optimal 
stopping problem similar to that for an American call option, but is subject to forced early exercise 
due to sudden job termination. From (12.81) . we can also interpret the vested ESO as an American call 
option with a cash flow stream of X(Soe Xt — K) + up to the exercise time r A . Using the ESO payoff 
structure, it is straightforward to show that the vested ESO cost C(t,x ) is increasing and convex in 
x for every t G [t v ,T), and is decreasing in t for every x € R. 

During the vesting period, the ESO is forfeited if the employee leaves the firm. Hence, given the 
ESO is still alive at time t <t v , the value of an unvested ESO is 

C(t,x) = E^| e~ r ^C(t v ,X tv )l {T - x>tv} | (2.9) 

= Eg J e~ ( r ' +A ) ( tv C(t v , X tv ) j. (2.10) 

The vesting provision prohibits the employee from exercising the option even if the ESO happens to 
be in the money during [0 ,t v ). 

The valuation of a vested ESO leads to the analytical and numerical studies of an inhomogeneous 
partial integro-differential variational inequality (PIDVI). To this end, we first define the infinitesimal 
generator of X under Q 

2 n 

Cf{x) = ilf'{x) + ^—f"{x)+ / (f(x + y) - f(x) - yf'{x)l {M<1} )i>(dy), (2.11) 

1 J R\{0} 

with fb given in (12.61) . For the vested ESO cost, job termination risk gives rise to an inhomogeneous 
term in the PIDVI, namely, 

mm^-(d t + C)C + {r + \)C-\(S 0 e x -K) + , C{t,x) - (S 0 e x - K) + j= 0, (2.12) 

for (t,x) G ( t Vl T ) x R, with terminal condition C(T,x) = ( Soe x — K) + , for id, 

For the unvested ESO cost, we set the terminal condition at time t v by matching it with the 
vested ESO cost, namely, C(t v ,x ) = C(t v ,x), for x € R. During the vesting period, the unvested 
ESO cost satisfies the partial integro-differential equation (PIDE) 

(dt + C)C - (r + X)C = 0, (f, x) G [0, t v ) x R. (2.13) 

When the vesting period coincides with maturity, the ESO becomes European-style as no early 
exercise is permitted. Setting t v = T yields the European ESO cost 

C E (t,x) = Eg, j e-^-^Soe^ - ^) + 1 {t a>t } }> 0 < t < T. 

This cost function also satisfies the PIDE (12.131) . 


>y 


J\-( r+x)( - u -VX{S 0 e x " - K) + du\, (2.8) 
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2.3 Exercise Boundary 

For a vested ESO, the holder’s exercise strategy can be described by the optimal exercise boundary 
t i->- s*(t) that divides the domain [t v ,T) x R into the continuation region C and exercise region S, 
defined by 


C = {( t,s ) € [t v ,T) x R + ,s < s*(f)}, ( 2 -14) 

5 = {(t,s) € [ t v ,T ) x R+,s > s*(t)}, (2.15) 

where 

s*(t) := sup{s > 0 | C(t, x) > (s — K ) + }, for t G [t v , T). (2-16) 

If s*(t) < +oo, then we have C(t,x) > (s — K) + for s < s*(t ), and C(t,x ) = (s — K) + for 

s > s*(t), for s = Soe x € (s*(f),+oo), due to the convexity and positivity of C(t,x ). Since, for every 
fixed x € M, C[t , x) is decreasing in t, the optimal exercise boundary s*(t) must be decreasing in t in 
view of ()2.16p . As for the impact of job termination risk, we have the following result: 


Proposition 2.1 A higher post-vesting job termination intensity decreases the costs of vested and 
unvested ESOs, and lowers the optimal exercise boundary. 


We pr ovide a proof in the Appendix IA.1I A similar result has been established under the GBM 
model by Leung and Sircar (2009). We remark that the cost reduction effect of the post-vesting job 
termination intensity holds with and without dividends. However, with q = 0, the job termination 
does not affect the employee’s voluntary exercise timing since it is optimal not to exercise early 
voluntarily regardless of job termination risk. In this case, the value of an American ESO equals to 
that of a European ESO, and the variational inequality (12.12|1 is simplified to a PIDE. On the other 
hand, a higher pre-vesting job termination intensity can reduce the unvested ESO cost, but has no 
impact on the vested ESO value or the post-vesting exercise strategy. 


Remark 2.2 In our paper, the optimal exercise boundary is computed based on the risk-neutral 
pricing measure. In practice, the ESO holder is likely unable to perfectly hedge the ESO exposure. In 
addition, the ESO holder may opt to exercise early due to other exogenous factors, such as liquidity 
risk or need for diversification. To this end, one can develop a reduced form or intensity-based 
approach by treating the ESO exercise timing as fully exogeneous, and calibrating to observed exercise 
behaviors. 


3 Fourier Transform Method for ESO Valuation 

For ESO valuation, we now discuss a numerical approach for solving the PIDVI (12.1211 based on 
Fourier transform. We first state the definition and some basic properties of Fourier transform. For 
any function f(x), the associated Fourier transform is defined by 

/ OO 

f{x)e~ wx dx , (3.1) 

-OO 

with angular frequency u in radians per second. In turn, if we denote by /(w) the Fourier transform 
of fix), then its inverse Fourier transform is 

~ 1 r°° ~ 

r 1[ / ] W = - J f(u)e™d». (3.2) 


8 

















As is well known, the Fourier transform of derivatives satisfies 


mf\H = iwjwr'fiw = ■■■ = ( 3 . 3 ) 

Applying (13.1 p to (12. lip , we have 

F[Cf]{u) = {ifiuj - \a 2 ur + [ (e^ M - 1 - iyul { \ y \ <1} )i'(dy)}J r [f}(uj) 

A JR 

= % >)Flf](u), (3.4) 


where 'F is the characteristic exponent under Q. 

In the continuation region, the vested ESO cost C(t,x ) satisfies the inhomogeneous PIDE 


-(dt + C)C + (r + X)C - tp(x) = 0, 


(3.5) 


where (p(x) : = X(Soe x — K) + . 

An application of Fourier transform to (13.51) yields 

d t T[C](t,uj) + (4f(u>) -r- \)F[C\(t,u) = -JM(u). (3.6) 


Therefore, the original inhomogeneous PIDE is transformed into an inhomogeneous ODE (13.61) sat¬ 
isfied by F[C]{t,uj), a function of time t parametrized by ui. Given the value of F[C\ at any time 
t 2 < T, we have at an earlier time t± that 

F[C](t i,w) = ^[C](f 2 ,o;)e^ (aj) - r - A)(t2 - tl) + ( . )( e (*M-r-\)(t 2 -*i) _ ^ ^ < ^ ( 37 ) 

’F(w) - r - A 


By inverse Fourier transform, we recover the vested ESO cost in the continuation region 

C(ti,x) = ^ 1 [^[C](t 2 ,u)e^^ ) - r - m2 - tl) }(x) + T~ 1 [{ , 

T(w) - r - A 

(3.8) 

Since the ESO is early exercisable, we need to compare the vested ESO value with the payoff from 
immediate exercise. Precisely, we partition the time interval [t v ,T] into t m ,m = M,M — 1,..., 1, 
then we iterate backward in time with 


C(t m - 1,X) = + 

C(t m -i,x) = max{C'(t m _i, x), (S 0 e x - K) + }, 


JF[tp\(u) r (jl(u>)—r—X)(tni—trr,.-T) 


!)](*), 


(3.9) 


where tM = T and to =t v . 

For numerical implementation, we discretize the original domain 0 = [t v , T] x M into a finite grid: 
{(t m ,x n ) : m = 0,1,... , M, n = 0,1,... , N — 1}, where t m = t v + mAt, and x n = x m i n + nAx, with 
A t = (T — t v )/M and Ax = ( x max — x m i n )/(N — 1). As most ESOs are granted at the money, it is 
natural to set the upper/lower price bounds x m j n = —x max equidistant from zero. With x max , M, 
and N fixed, we apply the Nyquist critical frequency u max = ir/Ax and set Aw = 2oj max /N. 

The continuous Fourier transform is approximated by the discrete Fourier transform (DFT) 

N—l N—l 

F[C](tmiW n ) « ^ C(t m ,x k )e~ lUnXk Ax = a n ^ C(t m , x fe )e^, (3.10) 

k= o fc=0 
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with a n = e~ lUnXrnin Ax. In (13.101) . we evaluate the sum J2k=o x fc) e_ ^ using the Fast Fourier 

Transform (FFT) algorithm. The corresponding Fourier inversion is conducted by inverse FFT, 
yielding the vested ESO cost C(t m ,x n ). Note that the coefficient a n will be cancelled in the process. 
After computing the vested ESO values, the unvested ESO cost is given by 

C{t,x) =JF- 1 [JT[(7](:r), for t<t v , (3.11) 

where C(t v ,x ) = C(t v ,x). Again, the Fourier and its inversion are implemented via FFT. 

We now provide some numerical results to illustrate the application of the pricing methods dis¬ 
cussed above. In Table [3] we summarize the ESO costs for different vesting periods t v = 0,2,4. As 
a numerical check, we compare with the ESO costs computed from a finite difference method (see 
Appendix IA.3p . and observe that the two numerical methods return very close values. 


Model 

t v 

= 0 

t v 

= 2 

t v 

= 4 


FST 

FDM 

FST 

FDM 

FST 

FDM 

GBM 

1.3736 

1.3730 

1.3822 

1.3816 

1.2365 

1.2360 

Merton 

1.4820 

1.4803 

1.4899 

1.4887 

1.3313 

1.3306 

Kou 

1.4566 

1.4558 

1.4648 

1.4646 

1.3091 

1.3104 

VG 

1.5584 

1.5595 

1.5816 

1.5811 

1.4131 

1.4139 

CGMY 

1.8409 

1.8411 

1.8532 

1.8535 

1.6484 

1.6490 


Table 3: ESO cost comparison. For each vesting period (t v = 0,2,4), the three columns contain ESO costs 
computed, respectively, by FSTC, FSTG, and the finite difference method. Common parameters: S$ = 
K = 10, r = 0.05, q = 0.04, A = 0.2, A = 0.1, T = 8. In the GBM model, a = 0.2. In the Merton model, 
a = 0.2, a = 3,/i = 0.02, a = 0.045. In the Kou model, p = 0.5, a = 0.2, a = 3 ,g+ = 50, rj- = 25. In the 
Variance Gamma model, (i = —0.22, a = 0.2, n = 0.5. In the CGMY model, C = 1.1, G = 10, M = 10, Y = 0.6. 
The grid sizes for FST methods are the same, with x max = 6, M = 2048, N = 32768. 


Figure [2] illustrates how the optimal exercise boundary changes with respect to job termination 
intensity A and stock price jump intensity a. As A increases from 0.1 to 0.3, the optimal exercise 
boundary is lowered. As the stock price jump intensity increases, the optimal exercise boundary 
moves upward. Also, we remark that in all cases the exercise boundary is decreasing in time t. 

The cost impact of job termination intensity and vesting period is demonstrated in Figure [3j 
As suggested by Proposition 12.11 a higher post-vesting job termination intensity reduces the ESO 
cost. Also, when the post-vesting and pre-vesting job termination rates are the same, the ESO 
cost decreases as vesting period lengthens (see Figure [3] (left)). However, if the post-vesting job 
termination rate (A = 0.2) is higher than the pre-vesting rate (A = 0.1), then it is possible that the 
ESO can first increase with vesting period. 


Remark 3.1 The model proposed by Cvitanic et al. ( 200A ) assumes that the ESO holder will volun¬ 
tarily exercise as soon as the log-normal stock price reaches an upper exogenous barrier. In addition, 
they also incorporate a vesting period and constant job termination intensity. Our current framework 
can also be adapted to their model. Precisely, one can numerically solve the PIDE 


{-dt + £)C + (r + A )C - A(S’ 0 e x - K) + = 0, 


(3.12) 


with the modified boundary condition C(t,x ) — C(t,L(t )) = 0, for (t,x) € ( t v ,T ) x (log(L(t)),+oo), 
and terminal condition C(T,x) = (Soe x — K) + . Here, the function L(-) is a given exponentially 
decaying barrier. An ESO cost comparison is provided in Tabled below. 
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Figure 2: (Left) ESO exercise boundaries for different job termination rates. (Right) ESO exercise boundaries 


for different stock price jump intensities. This example is based on the Kou model with parameters: Sq = K = 
10, r = 0.05 ,p = 0.5, cr = 0.2, q = 0.04, r/ + = 50, r?_ = 25, A = 0.1, T = 8. 




Figure 3: (Left) The ESO cost decreases as vesting period t v lengthens or as post-vesting job termination rate 
A increases. (Right) With A = 0.2 and A = 0.1, longer vesting can increase the ESO cost. This example is based 
on the Kou model with parameters: So = K = 10, r = 0.05, p = 0.5, a = 0.2, q = 0.04, q + = 50, p- = 25, T = 8 . 



Barrier ESO 

9 = 0 

European ESO 

American ESO 

Barrier ESO 

9 = 0.04 
European ESO 

American ESO 

L=125 

22.7792 

37.5435 

37.5435 

15.4209 

16.5753 

18.2484 

O 

LO 

r—1 

II 

26.8375 

37.5435 

37.5435 

17.4808 

16.5753 

18.2484 

L=9999 

37.5450 

37.5435 

37.5435 

16.5751 

16.5753 

18.2484 


Table 4: ESO cost comparison. Under the GBM model, the Barrier ESO ass umes exercise at an e xogenous 
decaying boundary, L(t) = Le at for t v <t <T with a = —0.02, as presented in ( Cvitanic et all 2008 . Table 2). 
In the zero dividend (9 = 0) case, the European and American ESO costs coincide, and they dominate Barrier 
ESO cost. Other common parameters: So = K = 100, r = 0.05, a = 0.2, A = 0.04, A = 0.04, t v = 3, T = 10. 


In Table |4j the barrier ESO corresponds to Case D in Cvitanic et al. ( 20081 ). i.e. the optimal 
exercise boundary is an expo nentially decay in g cur ve: L(t) = Le at , where a < 0. Accordingly, we 
can see: (1) The ESO cost in Cvitanic et al. ( 20081 ) is always underestimated if the ESO holder is 
allowed to exercise the ESO after the vesting period, because the exogen ous exercise b o undary is not 
generally the real optimal exercise boundary for the ESO holder. (2) In Cvitanic et al. ( 20081 ). when 
L becomes larger and larger, indicating that the ESO holder is unlikely to exercise the ESO according 
to exogenous exercise boundary, the ESO cost gets closer to the European ESO cost in our paper. 
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This case is also discussed in Carr and Linetskv ( 200C)h . (3) When q = 0, the value of European ESO 
is equal to the value of American ESO, and when q > 0, the value of European ESO is less than the 
value of American ESO. In general, our algorithm can also be modified to adapt other appropriate 
payoff at job termination. In the special case with zero payoff at job termination, the ESO can be 
interpreted as being forfeited at the time of departure. Alternatively, this can be considered as an 
American call option with default risk and zero recovery. 


4 Stochastic Job Termination Intensity 

As an extension to our ESO valuation model, one can randomize the job termination rate by defining 


T A = inf{t>0: f A (s,X s )ds > E], 

Jo 


where A (t,x) is a smooth positive deterministic function and E ~ exp(l) _L F A . This Markovian 
job termination intensity allows for dependence between r A and X w hile preserving trac t ability 
by not increasing the dim ension of the PIDVI. In related ESO studies, ICarr and Linetskvl ( 2000 1 
and Cvitanic et al. ( 20081 ) also consider this Markovian intensity approach to model job termina¬ 
tion an d exogenous exercis e, though they do not incorporate optimal voluntary exercise. In par¬ 
ticular, Carr and Linetskv ( 2000 ) consider the job termination intensity of the form: A (t,Xt) = 
A f + A e l| 5oe x i>A 'j. The second term is to model the early exercise due to the holder’s exogenous 
desire for liquidity, and it is constant if the ESO is in-the-money and zero otherwise. 

In our paper, we assume the job termination intensity functions before and after the vesting 
period to be affine in log-price, denoted respectively, by A(t, x) = ax + b and A (t, x) = ax + b, for some 
constants a, b, a, and b. Intuitively, the ESO holder’s employment is more at risk when the company 
stock price is low, so one can let a and a be negative. For implementation, one can select parameters 
and control the grid size so that the intensity remains positive within the truncated log-price interval 
The affine intensity assumption is utilized in simplifying the associated inhomogeneous 


[■f 'm.in j %max 

PIDE. 


The PIDE for the vested ESO cost in the continuation region is 

(dt + C)C — rC = (ax + b)C + ip(x), 


(4.1) 


where 'ip(x ) = (ax + b)(Soe x — K) + . Differentiating (jf.ip w.r.t. x, and applying Fourier transform 
and (|3.3D . we obtain an inhomogeneous PDE in t and cj, namely, 


iui(dt + ('k(w) — r)F[C](t,uj)) = F[aC + (ax + b)C x \(t,uj) — iui F[ip\ (oj) , 

with terminal condition F[C](T, oj) = J r [(5oe x — K) + ], 

Using the following well-known property of Fourier transform: 

F[xC x \(t,u) = -T[C](t,u) - ud„F[C](t,u), 

followed by the substitution 

E[C\(t,u) = e^f-~^- r - b)ds H{t,uj), 
we simplify (14.21) to a first-order PDE 

d t H(t,uj) + -d u H(t,uj) = 


(4.2) 

(4.3) 

(4.4) 

(4.5) 
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with terminal condition H(T,uj ) = F[(Soe x — K) + )(uj)e a f-oo(^(0 1 b ) ds . Finally, solving (14. 5p and 
unraveling the substitution, the Fourier transform of C can be expressed as 




u ,_ l _ A /* A£ 

J-[C(ti,x)](w) = e^^-- (vl>(s) - r - b),is (^[C(t 2 ,x)](u;)e^—“ (^( s )-r—6)rf s + / “ ( 5 _ At))ds ), 

JO 1 

(4.6) 


where 


<?(w) = e a f-oo(^( s ) r b ^ ds F[f>\(u), t v <t\<t 2 <T, At = t 2 — t\. 


The Fourier transform (14. 6p allows us to compute the values of C(t, x ) backward in time, starting 
from expiration date T. The numerical implementation of (|4.6I) requires computing the integral 
/ 0 A * g{ui — f(s — At))ds. Within each small time step [t, t+ At], we approximate the integral J 0 A< g(ui — 
j(s— At))ds by summing over a further divided time discretization, namely, 3 ^ F At) At / n'. 

Therefore, the solution for vested ESO cost in the continuation region is given by 


C(h,x) 


n'—l 


= F~ l [e- ^ At ^- r - b)ds T[C{t 2 ,x)\{u + -At) + y g(u + ^-At)^-e^ ( x ), 

L i m n 

k =o 


for t v < t\ < t 2 < T, x € R. Again, within each iteration, we impose the condition C(t\,x) = 
max{(7(ti, x), ( Soe x — K) + } to get the American option value. 

As for the unvested ESO, we have C(t v ,x ) = C(t v ,x ) at time t v , and 

(dt + C)C - (r + (ax + b))C = 0, 


for (t,x) € [0,t„) x R. At time t < t v , the unvested ESO cost is computed by 


C(t,x ) = F 1 [. F[C(t v ,x)\(uj + ~(t v — t))e a ^ 




(4'(s)-r-b)dsj 


(4.7) 


Remark 4.1 Suppose the pre-vesting and post-vesting job termination intensities, A and \, are pos¬ 
itive bounded functions of x and t. The Fourier transform of the ESO cost satisfies 

dtT[C](t,u) + (tf (w) - r)F[C](t,u) - F[XC](t,co) = -.F[£]M, (4.8) 


where £(t,x) = X(t,x)(Soe x — K) + . In order to solve this ODE, one can apply an explicit scheme 
to the term F[XC](t,oj), and apply implicit scheme to the other terms. Therefore, given the value of 
C(t 2 ,x) at time t 2 > t\, we compute the value of C(t\,x) by 


C(t!,x) = F- l [E[C\(t 2 ,w)e ( ’^ ) - r)(t2 - tl) ] + F~ x jr i XC \( t ^ u ) 

T(u;) — r 


With this, we can compute the vested ESO cost by iterating backward in time and comparing with 
the payoff from immediate exercise. Finally, we remark that this implicit-explicit algorithm is also 
applicable when the job termination intensity is constant or affine. 


In Table 0 we compute the ESO cost with affine pre-vesting and post-vesting job termination 
rates using the FST method in Section [4] and that in Remark 14.11 along with the finite difference 
method in Appendix IA.3I We observe that ESO costs for different vesting periods from these three 
different methods are very close. 
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Model 


t v = 0 



t v = 2 



tv = 4 



FSTA 

FSTG 

FDM 

FSTA 

FSTG 

FDM 

FSTA 

FSTG 

FDM 

GBM 

1.3732 

1.3733 

1.3729 

1.1368 

1.1369 

1.1364 

0.8429 

0.8428 

0.8429 

Merton 

1.4817 

1.4814 

1.4800 

1.2261 

1.2259 

1.2260 

0.9085 

0.9083 

0.9082 

Kou 

1.4565 

1.4562 

1.4556 

1.2054 

1.2052 

1.2064 

0.8934 

0.8933 

0.8942 

VG 

1.5670 

1.5669 

1.5680 

1.3073 

1.3069 

1.3068 

0.9765 

0.9754 

0.9758 

CGMY 

1.8402 

1.8398 

1.8402 

1.5249 

1.5247 

1.5245 

1.1299 

1.1297 

1.1296 


Table 5: ESO cost comparison under affine job termination intensity. For each vesting period (t v = 0, 2,4), the 
three columns contain ESO costs computed, respectively, by the two FST methods in Section [4] and Remark 
14. H and a finite difference method. The job termination intensities are \{x) = A(x) = —0.022 + 0.2, while other 
parameters are the same as in Table [3] 

5 Risk Analysis for ESOs 

Over the life of an ESO, the option value will fluctuate depending on the company stock price 
movement. From the firm’s perspective, an increase in the ESO value implies a higher expected cost 
of compensation as compared to the initially reported value. For the purposes of risk management 
and financial reporting, it is important to consider the probability that the ESO cost will exceed a 
given level. 

In addition, the ESO can be terminated early voluntarily or involuntarily prior to the expiration 
date. This also motivates the study of the contract termination probability, which can be a tool to 
calibrate the job termination intensity parameter. Such a calibration would provide a crucial input 
to the ESO valuation model, and permit the pricing to be consistent with the firm’s characteris¬ 
tic. Moreover, it is also interesting to examine the impact of job termination risk on the contract 
termination probability. 

Generally speaking, the job termination rate can differ under P and Q. However, since ESOs are 
not traded, market prices are unavailable for inferring the Q intensity. For this reason and notational 
simplicity, we assume that the historical and risk-neutral job termination rates are identical. 

5.1 Cost Exceedance Probability 

Let us evaluate at the current time t the probability that the ESO cost will exceed a given level £ on 
a fixed future date T . To this end, we need to consider three scenarios separately. 

Case 1: t < T <t v 

In this case, the time interval [t,T] lies within the vesting period [0,i„], where the unvested ESO 
can be forfeited upon job termination. Therefore, the probability that the ESO value will exceed 
a pre-specihed level £ > 0 at time T is given by P t>x {C(T,Xf) > £}, where Pt jX is the historical 
probability measure with Xt = x. Since the ESO becomes worthless if t x arrives before T, we have 

p t,x{C(T, x f ) > £} = P t, x {C(T, Xf) > £, r x > T} + P t, x {C(T, Xf)>£,r x <T} (5.1) 

= P t , x {C(T, X f )>£\ r x > T}F ttX {r x > T}. (5.2) 

In (15.21) . we notice that P t,x{T X > T} = e~ x ^ T ~^ since t x ~ exp(A). Hence, the evaluation of 
P t,x{C(T 1 Xf) > £} amounts to computing the conditional probability P t lX {C(T, Xf) > £ \ r x > T}. 
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Since C(T,x) is increasing in x, we can find the critical log-price x such that C(T,x) = £ and 
write 

p(t, x) := P t, x {C(T, X f )>£\ r x >T} = P tiX {X f >x\t x >T}. (5.3) 

The probability function p(t, x ) satisfies the PIDE problem 

(dt + C)p = 0, (t, x) € [0, T] x M, (5.4) 

p{f,x) = l {x>s} , xeM, (5.5) 

where C, is the infinitesimal generator of X under P. We observe that the PIDE (15.41) for p(t, x) is 
very similar to (13.51) without the inhomogeneous term. Applying the Fourier transform arguments 
from Section [3] yields that 


p(t,x) = { , >s} ](w)e^^( f - t )](x) J t <T <t v . (5.6) 


The Fourier transform and its inversion in (15.61) can be numerically evaluated by the FFT algorithm. 
In contrast to the ESO cost, this probability does not involve early exercise and can be computed in 
one time step. 


Remark 5.1 An equivalent way to obtain (15.61) is to adapt the convolution method used by Lord_ et a,l 
(200 A) . To verify this, we express «op as 


P(t,x) 


l{z>x}fx f \x t (z)dz 



^{x-\-y>x} fXj,_ t {y)dy, 


(5.7) 


where fx \x t ( z ) °f Xf is the conditional probability distribution function given X t = x. Then, we 
express (15.71) in terms of Fourier transform: 


/ OO 

-oo 



/ +oo 

1 {x+y>x}fx t _ t (y)dy)dx 
-oo 

/ +oo 

1 {u>x}fx f _ t (y)dy)du (u = x + y) 

-oo 

/ +oo 

e iuy fx t _ t {y)dy = .F[l {a;>s} ]Me*H(T-t). 

-oo 


Lastly, applying inverse Fourier transform to the last equation yields ra, and hence the equivalence. 


A slightly different approach is to adapt the Fourier transform method by Carr and Madan ( 19991 ). 
To illustrate, we assume without loss of generality that Xt = 0 at a fixed time t < T, and define 


P(t,z) := Pt,o{X f > z\t x > T} = / fx f (y)dy, 


s: 


(5.8) 


where fx, r -_ t is the probability density function of Xf_ t with X$ = 0. The value z can be considered 
as the difference between the upper threshold x and the current value of X. 

Applying Fourier transform to both sides of (15.81 ). we get 

/ oo roo poo ry 

e~ ZU}Z ( f X j,_ t (y)dy)dz = / fx f _ t {y){ e~ luz dz)dy. 

-oo J z J —oo J —oo 
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We notice that the inner integral J^ oo e~ luz dz = e _^~ \ y _^ does not converge, so we incorporate a 
damping factor e az , for a > 0, and consider p a (t,z) = e az p(t,z). Then, the corresponding Fourier 
transform is given by 

/ oo ry 

fx t _ t (y ) e {a ~ luj)z dzdy = (5.9) 

-oo J—oo a ^ 

where is the characteristic function of Xf_ t . In turn, inverse Fourier transform to (15.91) yields 


P(t, z) 


e az 

2ir 



e lzU} <l)f_ t (—uj 
a — ioj 




e luz ,— ia ) 
a — iuj 


du, 


(5.10) 


where we have used the fact that the Fourier transform of the real function p(t, z) is even in its 
real part and odd in its imaginary part. Lastly, the integral in (15.101) is approximated by the FFT 
algorithm. We remark that p(t,z) is independent of the choice of a, and we choose a S [1.5, 2.0] for 
numerical implementation. 

In Table [6j we present the numerical results for the ESO cost exceedance probabilities under dif¬ 
ferent thresholds l. The FST-FST method first computes the ESO cost vector via FST (see (13. lip ) 
which gives the critical value x fsee 15.31) . In the second step, it applies FST to solve for the exceedance 
probability as shown in (15.61) . The FFT-FST method differs from the FST-FST method in its second 
step where FFT (see (15.101) 1 is used to compute the probability. For the GBM, Merton, and Kou 
models, the reference values are given by closed-form formulas (see Appendix B.2). The numerical 
results show that both Fourier transform methods are very accurate. 


Model 

1 = 

Reference Value 

1.1 C(0,0) 
FST-FST 

FFT-FST 

£ = 

Reference Value 

1.2 C(0,0) 
FST-FST 

FFT-FST 

GBM 

0.2534 

0.2532 

0.2535 

0.0868 

0.0869 

0.0869 

Merton 

0.2607 

0.2608 

0.2607 

0.0943 

0.0942 

0.0943 

Kou 

0.2431 

0.2430 

0.2431 

0.0799 

0.0797 

0.0798 


Table 6: Cost exceedance probability in case 1. Here, t = 0,T = 10/252, t v = 2,T = 8 and other parameters 
are the same as in Table [3] 


Case 2: t v < t < T < T 

After the vesting period, the employee intends to exercise at any t* < T, but may be forced to 
exercise at r A . The probability of interest is 

p(t, x ) := Pt,s{C(f* A r\ X f * Ar x) > £}, (5.11) 

with r* = t* AT. To compute this, we first identify the critical log-price x{t), such that C(t, x{t)) = £, 
for t € [0, T], In turn, we write down the corresponding inhomogeneous PIDE 

(<9i + C)p — Xp + A1{ x > S (£)} = 0, (5-12) 

for (t, x) € ( t v ,T ) x M. The boundary and terminal conditions depend on the relative positions of the 
critical log-price x(t) and optimal exercise boundary x*(t) := log(s*(t)/S'o). Precisely, if x(t) < x*(t), 
then we set p(t,x) = 1 at time t for x > x*(t). If x(t) > x*(t), we set p(t,x ) = 1 for x > x{t), and 
p(t,x ) = 0 for x*(t) < x(t). As for the terminal condition, if x(T) < K, then we set p{T,x ) = 1, for 
x > K. If x(T) > K, then we have p{t,x ) = 1 for x > x(t), and p(t,x ) = 0 for K < x < x(t). 
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For numerical implementation, we first solve for the ESO values C(t,x), which are then used to 
determine the critical log-price x(t), and the associated optimal exercise boundary s*(t) that gives 
x*(t). With these, we apply Fourier transform implicit-explicit method, as discussed in Section [4.11 
to the PIDE problem (15.121) . Iterating backward in time, we have for each time step 


p(t m - i,x) 

v (cj j — A 


!)](*), ( 5 - 13 ) 


for m = M ,..., 1,0, along with the boundary conditions. In Figure [4] the cost exceedance probability 
rises as the horizon T lengthens. The probability also increases as the job termination intensity A 
decreases. 



Figure 4: ESO cost exceedance probability in case 2 increases with T, and decreases with job termination rate 
A. This example is based on the Kou model with common parameters: S 4 = K = 10, r = 0.05, p = 0.5, q = 
0.04, p = 0.08, A = 0.1, T = 8 , t = 4, a = 0.2, a = 3, = 50, r]_ = 25 


Case 3: t < t v < T < T 

This scenario is a combination of case 1 and case 2. The ESO is forfeited if r A < t v , as in case 1. 
However, if that does not happen, then case 2 applies after the vesting. Consequently, we consider 


the following probability for ESO cost exceedance 

p(t,x) := EjTjl {T% > tv} p(t v ,X tv )}, 0 <t<t v , (5.14) 

if the job termination has not occurred by time t. The corresponding PIDE problem is 

(d t + C)p — Xp = 0, (t, x) G [0, t v ) x R, 

p(t v ,x) = p(t v ,x), x G R. (5.15) 

The solution via Fourier transform is given by 

p(t v ,x) = p(t v ,x), (5.16) 

p(t,x) = J 7 ~ 1 [j 7 \p(t v ,x)](cj)e^^~^^ tv ^](x), for t<t v . (5-17) 


Here, the probability p(t v ,x ) from case 2 is used as the input, and then p(t,x ) can be computed 
directly without recursion by (15.171) for any time t <t v . 
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5.2 Contract Termination Probability 

The ESO contract can be terminated either due to job termination before vesting, or voluntary/forced 
exercise after vesting. We now study the ESO contract termination probability during a given time 
interval [t, T]. Again, we consider three different scenarios. 

Case 1: t <T <t v 

During the vesting period, the contract termination is totally due to job termination before t v , of 
which the probability is simply 1 — e _A ( T ~^. 


Case 2: t v <t < T <T 

After vesting, contract termination can arise from either involuntary exercise due to job termina¬ 
tion, or the holder’s voluntary exercise. For calculation purpose, we divide the contract termination 
probability into two parts according to whether job termination occurs before or after T. 

First, we consider the scenario where job termination does not occur during [t, T] and the ESO 
holder voluntarily exercises the ESO. This corresponds to the probability 

P t, x { T * <T,t x >T} = P \ x {t* < T}e- A M, (5.18) 

where r* is the holder’s optimal exercise time. To evaluate (15.181) . we solve for h(t,x ) := P t,x{T* < T} 
from the PIDE 


(dt + C)h — 0, 


(5.19) 


for (t, x ) € ( t v , T) x M, with the boundary condition h[t , x) = 1, for x > x*(t), and terminal condition 
h(T,x ) = where x*(t ) := log(s*(t)/5o). For numerical solution, we apply recursively 

[h(t m -i,x) = J r_1 [J'[/i(t m ,a:)](a;)e 5 ' (tj)(<m “ t " 1 - l) ](x), x<x*(t) 

lh(t m -!,x) = 1, x > x*(t). 

The other scenario is when job termination arrives before T. Consequently, the total contract 
termination probability h(t,x ) is the sum 

h(t,x) :=e- A ( T_t) ^(t,s) + l-e- A ( T_t ). (5.21) 

From this expression, we observe that the contract termination probability is increasing with job 
termination intensity A since the optimal exercise boundary is decreasing with A, and so is (1 — h(t, x)). 
Alternatively, we look at the employee’s voluntary exercise probability 

h v (t,x ) := Pt,x{r* < r A A T}. 

This probability satisfies the PIDE 

(d t + C)h v - A h v = 0, (5.22) 


for (t, x) € ( t v ,T) xR, with the boundary condition h v (t, x) = 1, for x > x*(t), and terminal condition 
h v (T,x) = 1 { a .> a ,*(f’)}, where x*(t) := log(s*(t)/S q). The numerical solution is found from backward 
iteration with 


h v (t m -i,x) = J'- 1 [J r [h ,; (t m ,x)](a;)e( vl> ^)- A )( t -- t "‘- 1 )](x), x < x*(t) 
h v (t m -i,x) = 1, x>x*(t), 
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for m = M, ...,1,0. 

Figure [5] illustrates the contract termination probability on top of the voluntary exercise proba¬ 
bility. From PIDE (15.22|) . we observe two competing factors governing the effect of job termination 
intensity A on the voluntary exercise probability. On one hand, a higher job termination intensity 
A implies a lower optimal exercise boundary, which in turn increases the voluntary exercise prob¬ 
ability. On the other hand, a higher job termination intensity will more likely force early exercise 
before the stock price reaching the holder’s optimal exercise boundary. This reduces the voluntary 
exercise probability. In Figure [5] (left), we see that voluntary exercise probability is decreasing with 
post-vesting job termination intensity A, so in this case the job termination effect outweighs the effect 
of lowered exercise boundary, and therefore, reduces the probability of voluntary exercise. Also, in 
Figure [5] (right), the contract termination probability is increasing with job termination intensity A, 
as expected. 




Figure 5: (Left) Job termination intensity increases contract termination probability, but decreases voluntary 
exercise probability. (Right) Stock price increases both contract termination probability and voluntary exercise 
probability with A = 0.1 and A = 0.2. This example is based on the GBM model with common parameters: 
S 6 = 10, K = 9, r = 0.05, cr = 0.2, q = 0.04, n = 0.08, t v = 2, T = 8, t = 6, f = 7. 


Case 3: t <t v < T <T 

This scenario is a combination of cases 1 and 2 above. The contract termination can occur before or 
after vesting. During [t, t v ], only job termination can cancel the contract. If there is no job termination 
before t v , then the contract termination resembles that in case 2. Therefore, the contract termination 
probability is the sum 


! _ e -A(G-t) +e p x {h(t v ,X tv )l {r ^ tv} } 


=:h(t,x) 


(5.24) 


where h(t,x ) is given in (15.2111 . Hence, h(t,x ) satisfies, for ( t,x ) € [0 ,t v ) x R, the PIDE 

{dt + C)h — Xh = 0, 


(5.25) 


At time t v , we set h(t v ,x ) = h(t v ,x), where h(t v ,x ) is computed from case 2. Again, the FST method 
discussed in Section [3] can be applied to solve for h(t,x). At any time t < t v , the probability can be 
computed in one step (without time iteration) via 

h(t,x) = J r_1 [J r [/i(t„,a:)](u;)e ( ^ (aj) " A)(t ’'" t) ](x). (5.26) 
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6 Perpetual ESO with Closed-Form Solution 


We now discuss the valuation of a perpetual ESO under the GBM model. In contrast to the valuation 
under the general Levy framework, the perpetual ESO admits a closed-form solution, and thus 
represents a highly tractable alternative. We first consider a vested perpetual ESO, whose value 
can be expressed in terms of an optimal stopping problem, namely, 


V(s) = sup 
reTo.oc 


e ~(r+X)T( Sr _ jq+ + J e -(r+A > X ( Su _ K )+ du | So = 4, 


where St is the real stock price. The associated variational inequality is 

mini -CV + (r + \)V - \{s - K)+, V{s) - (s - K) + 1=0, s € R+. 


( 6 . 1 ) 


( 6 . 2 ) 


A similar inhomogeneous variational inequality has been derived and solved for the problem 
of pricing American puts with maturity randmization (Canadization) introduced by Carr ( 19981 ). 
Indeed, the perpetual ESO can be considered as an American call whose maturity is an exponential 
random variable. For a mathematical analysis on t he C anadization of American options with a Levy 
underlying, we refer to Kvprianou and Pistoriusl ( 2003 ). Next, we present the closed-form solution 
for this ESO valuation problem. 

Proposition 6.1 Under the GBM model, the value of a vested perpetual ESO is given by 

' Dsi+ 


V{s) = 


As 1+ + Bs^~ + j^s - if K < s < s* 

s-K 


if s < K, 
if K < s 
if s > s* 


(6.3) 


where 


7± = 
B = 
A = 


(q - r + 4 ) ± \J(q - r + ^) 2 + 2(r + X)a 2 


A(1 - 7 +) 


a* 


7 _( 7+ - 7 _)(A + g) 

4 (s*) 1 ”^ - ^—B (s*) 7 — 7+ , 

(A + q) 7 + 7 + 


D = A + B + -44—- K l ~ 1+ . 

(A + g)(A + r) 

The optimal exercise threshold s* € (K, oo) is uniquely determined from 


B(1-— )(8*) 7 -- (1 - L)( 

7 + 7+ A + q 


Q )s* + ^-K = 0. 
r + A 


(6.4) 

(6.5) 

( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 


A number of interesting observations can be drawn from Theorem (16.11) . First, in the case without 
job termination (A = 0), we have B = 0 and A = D. This implies that the perpetual ESO reduces to 
an ordinary American call with the well-known price formula 


V(s) = 


j Ds' 1+ if s < s *, 
| s — K if s > s *, 


(6.9) 
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where s* = 1 2^ + K. This result dates back to Samuelson j 19651) and McKean ( 1965 1. and is also used 
in the real option literature (see, for example, McDonald and Siegel ( 1986l l). Furthermore, if both A 
and q are zero, then we can see that s* = oo, which means that it is optimal for the option holder to 
never exercise the option. This is expected since the ESO now resembles an ordinary American call 
without dividend. 



Figure 6: As job termination rate increases (A 6 {0,0.2,1}), the vested ESO cost decreases and the correspond¬ 
ing optimal exercise threshold is lowered (sj$ > Sq 2 > s i)- Parameters: K = 10, r = 0.05, a = 0.2, q = 0.04. 


When a vesting period of t v years is imposed, we compute the ESO cost at time t from the 
conditional expectation 


V(t, s ) = E Q {e -(r+A)(t » -t V(StJ | St = a}. 


( 6 . 10 ) 


By substituting the vested ESO cost function V(s) into (|6.10j) . and recognizing that , for any 
a G R, is lognormal, we can directly compute the unvested ESO cost. 


Corollary 6.2 Under the GBM model, the unvested perpetual ESO cost admits the formula 
V(t, s) = P W"+0+ Ml + d) 


( 6 . 11 ) 


+ As 7 +e( mi+ ir) 

cr 2 r 

+ Ss 7 -e (m2+ ^ ) 


7 + In (s/K) + m\ + (t\ 


0 1 


7- ln(s/s*) -)- m 2 + a\ 


02 


- <b 


- $ 


7 + ln(s/s*) + mi + af 


01 


7 _ \n(s/K) +m 2 + cr\ 


02 


+ »(1 + _ 2_) e <"»+3 >* ( Ws / K )+™ 3 + 4 \ _ (1 

A + q J V 03 J 


r + A 


)K& 


In (s/K) + m 3 + cr| 


03 
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where $ is the standard normal complementary c.d.f. and 


m i = (r- 

™2 = {r-—)(t v 
, a 2 

m3 = ( r-—){t v 


t)l+, cri = oy/t v -t 7+, 
t) 7-i cr 2 = -ay/t v -t 7_, 
t), (J3 f ■ 


( 6 . 12 ) 

(6.13) 

(6.14) 


In contrast to its vested counterpart, the unvested ESO cost is time-dependent. In Figure [7] (left), 
we observe that the perpetual ESO cost decreases as the vesting period increases. A higher job 
termination rate not only decreases the ESO cost (see Figure [6]) but also the exercise threshold s* 
(see Figured] (right). 




Figure 7: The perpetual ESO cost decreases as vesting period increases (left) and as job termination rate 
increases (right). Parameters: Sq = K = 10, r = 0.05, a = 0.2, q = 0.04, t = 0 (left), t v = 0 (right). 


7 Conclusions 

We have provided the analytical and numerical studies for the valuation and risk analysis of ESOs 
under Levy price dynamics. Our results are useful for reporting ESO cost, as mandated by regulators, 
and for understanding holder’s exercise behavior. In particular, we show job termination risk has 
a direct effect on the ESO holder’s exercise timing, which in turn affects the ESO cost as well as 
contract termination probability. For future research, risk estimation for large ESO portfolios is both 
practical and challenging. Other related issues include the incentive effect and optimal design of 
ESOs and other compensation schemes, such as restricted stocks. Lastly, our valuation framework 
can also be applied to pricing American options with liquidity, default, or other event risks. This 
would require an appropriate modification of the payoff at the exogenous termination time. 
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A Appendix 

A.l Proof of Proposition 12.ll 

Let Ci(t,x) and C 2 {t, x) be the vested ESO cost associated with Ai and A2, respectively, and assume 
Ai < A2. We define an operator M by 

MiC{t, x ) = [d t + t)C - (r + A i)C + Xi{S 0 e x - K)+. (A.l) 

From the variational inequality (|2.12)l . we see that MiCi < 0. We choose a point (t,x) in the 
continuation region of C 2 , which means that M. 2 C 2 = 0. Since Ai < A2 and C 2 > ( Soe x — K) + , direct 
substitution shows that M. \ 62 > 0. 

Next, we define the process 

m(t,X t ) = e~( r+Xl)t C 2 {t,X t ) + [ g-h+^Ai (S 0 e Xu - K) + du, t> 0. (A.2) 

Jo 

Using the fact that M.\C 2 > 0 and Optional Sampling Theorem, we deduce that, for any r € lt,T, 

E ?x{ m (uA T )} > m(t,x). (A.3) 

In particular, we denote rj* and r 2 as the optimal stopping times associated with C 1 and C 2 , and get 

C 2 (t,x) < E ^|e _(r+Al){r2 *“ t) C' 2 (r 2 *,X T |) + £ 2 e-^^^Ad^oe^ - iL) + dn| (A.4) 

= - A) + | (A.5) 

<Ci(t,x). (A. 6) 

The last inequality follows since r 2 is one candidate stopping time for the optimal stopping value 
function C\. Hence, we conclude that C\ (f, x) > C*2(t, x) > (Soe x — K) + . This implies that any point 
(t,x) in the continuation region of C 2 must also lies in the continuation region of C\, which means 
that the optimal exercise boundary for C\ dominates that for C*2. 

As for the unvested ESO, the job termination intensity reduces its terminal values, and increases 
the probability of forfeiture (with payoff zero) during the vesting period. As a result, a higher job 
termination intensity also reduces the unvested ESO cost. 

A.2 Proof of Proposition IHTlI 

We conjecture that it is optimal to exercise the ESO as soon as the stock reaches some level s* > K. 
Then, we split the stock price domain into three regions: [s*, 00), [K, s*), and [0, A). In region 1, we 
have s > s* and V(s) = s — K. In region 2, the ESO cost solves the inhomogeneous ODE 

2 2 

^-V'\s) + (r - q)sV'(s) - (r + A)U(s) + A(s - K) = 0. (A.7) 

One can check by substitution that the general solution to (1A.7I) is given by 

V(s) = As 7 + + Bs 1 - + — s - ^—rK, (A.8) 

A + q r + A 
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where 7_ and 7+ are given in (16.411 . 

In region 3, since the option is out of the money, we have the ODE 

2 2 

^-V"(s) + (r - q)sV'(s) - (r + A)V(s) = 0 , (A.9) 

whose solution is the form V(s ) = Ds 7 + + Es 1- . Since h(s) -> 0 as s -> 0, it follows that E = 0. 

To solve for the constants A, B, D, along with the critical stock price s*, we apply continuity and 
smooth-pasting conditions at s = K and s = s* to get 


lim V(s) = V(K) 

s^K 

=4- 

PV+ = AK^+ + BK^~ + (, A A , )K, 

A +q r+ y 

(A. 10 ) 

lim V\s) = V\K) 

s^K 

=4 

7 + pr +" 1 = 7 + AAr 7+_1 + 7 _sr - _1 + A , 

A + q 

(A.ll) 

lim V(s) = V(s*) 

sfs* 

=4 

A(s*y++B(s*y-+ x X s* X K = s* K, 

A +q r+ A 

(A. 12 ) 

lim V'(s) = V'(s*) 

sfs* 

=4 

7 + A(V ) 7 +- 1 + t -Bisy- 1 + = 1 . 

A q 

(A. 13) 


Solving this system of equations yields (16.5D - (16.611 . In particular, we see that B > 0 since 7 + > 1 and 
7 - < 0 . 

From (I A. 12 II - (1 A. 13 f) . the threshold s* satisfies the equation 

/(»•) := B( 1 - - (1 - 4-)(^_) s * + -J—K = 0. (A.14) 

7+ 7 + A + q r + A 

To show this has a unique real solution, we note that / is continuous, and 

/'(«*) = 7 -B(l - —)(s *) 7 - -1 + (1 - —)(_Z5_) < 0, (A.15) 

7+ 7+ A + q 

since B > 0 , 7 _ < 0, and 7 + > 1. In addition, we have the limits: lim s *|o = 777 > 0 and 
lim^-fvx, = — 00 , as well as f(K ) = > 0. Together, this implies f(s*) = 0 has a unique real root 

s* > K. Hence, we obtain formula (16.311 for V. By direct substitution, it satisfies the VI (J6.2H . 


A.3 Finite Difference Method for ESO Valuation 

We summarize the finite difference method (FDM) for computing the ESO costs i n Tables [3] and [H 
For th is purpose, we adapt the FDM algorithm for European options detailed in Cont and Voltchkova 
( 2003h to the current case of early exercisable ESO with job termination risk. 

First, we introduce the change of variable u = T — t and denote F(u, x) = C(T — t,x). Then, the 
PIDE for the vested ESO cost in the continuation region becomes 


BF 

- 7 — = CF — (r + A(T — u, x))F + A(T - u, x)(S 0 e x - K)+, 

Bu 

for (u,x) € (0,T — t v ) x R, with the initial condition F(0,x) = (Soe x — iv) + , x € R, where 


(A.16) 


~ , ,BF a 2 ,3 2 F 8 F x 

CF - {r - q) te + T ( a? _ & ) + 


v(dy)(F(u,x + y) - F(u,x) - (e y - l)4fj. (A. 17) 

Bx 
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To proceed, we split the operator £ into two parts, namely, 

CF = DF + JF, 


(A.18) 


where 

_2 o2 77 ^_2 ar rB T 

= y-^ 2 -- (y - r+ g + /3)—- aF, and JF = J v(dy)F(u,x + y), (A.19) 

with /3 = /^ r i>(dy)(e y - 1). 

We define a uniform grid on [0, T — t„] x [—A, A] by {(u n , Xi) : u n = nAt, n = 0,1, 2,... , M, Xi = 
— A + iAx,i € 0,1,2..., A^}, with At = (T — t v )/M, Ax = 2A/N. In Tables [3] and 0 we have 
M = 4000 and N = 8000. Also, denote F™ be the cost at the grid point ( u n ,Xi ). We use trapezoidal 
quadrature rule with Ax to approximate the integral terms in (1A.19I) . To do so, we let Ki , K r be 
such that [Bi,B r \ C [(Ki — l/2)Ax, (K r + l/2)Ax], and apply the approximations 


rBr Kr K r j ~ R . 

/ 0(dy)F(u,Xi + y) « ^ VjF i+j , a « ^ Oj, (3 « ^ 

j=Ki j=K t j=K l 


uAe 


Vj _ 


1 ), 


with 


/• (.7+1/2) Ax 


v(dy). 


The space derivatives are approximated by the finite differences 


<9F <9 2 F _ - 2q + Fj_ 1 

<9x * Ax ’ <9x 2 * (Ax) 2 


(A.20) 


Next, we replace DF and JF with their approximations D A F and J A F , respectively. Lastly, we 
arrive at the following implicit-explicit time-stepping scheme: 


pn +1 _ p 


n 


At 


where 


D A F n+1 + J A F n - (r + A(T - (n + 1) At, x))F n+1 + A(T - (n + 1) At, x) ( S 0 e x - K)+, 

(A.21) 


(D A F)i 


c 2 Fi + 1 — 2Fi + T)_i 

~2 (Ax) 2 


X , , a\ F i+l~ F i r 

( __ r+9+ffl __— aF „ 


K r 

(JA F )i — ^ ' VjFi + j. 
j=Ki 


(A.22) 


Due to the early exercise feature, the iteration is coupled with a comparison to the payoff from 
immediate exercise. After computing the vested ESO cost till the end of the vesting period, similar 
finite difference method can be applied to solve the PIDE for the unvested ESO cost: 


dF 

du 


CF — (r + A(T — u, x))F, 


(A.23) 


for (u,x) € (0 ,t v ) x R. 

The above algorithm works for the case when the underlying Levy process has finite activity, with 
z>(R) = a < + 00 . In the infinite activity case with £>( R) = + 00 , we we can use an auxiliary process 
0 with the Levy triplet (/i(e),cr 2 + a(e) 2 ,01\ x \ >£ ) to approximate the original process (Xt)t>o, 
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where cr(e) 2 = ft^y 2 k'(dy), and /t(e) is again determined by the risk-neutral condition. Therefore, 
F £ satisfies 

8F £ 

-r- = C E F £ - (r + A(T - u,x))F £ + A(T - u,®)(,Soe* - K)+, (A.24) 

au 

for (u,x) € (0, T — t v ) x R with initial condition F £ (0,x) = ( Soe x — K) + , igl. Here, the operator 
C £ is defined by 


t £ F £ = ( 


with 


cr 2 + cr(e) 2 d 2 F £ a 2 + cr(e) 2 
2 ' dx 2 ~ ^ 2 


dF £ 

-r + q + - a(s)F £ (x) + 


>£ 


v{dy)F e {x + y), 
(A.25) 


/3(e) = f (e y — l)z>(ch/) and a(e) = f v{dy). 

J\y\> E J\y\>£ 

The PIDE in (IA.24D can be solved by the same numerical scheme as in the finite activity case. 
We apply this finite difference method for comparing with our FST method. For alternative finite 
difference methods, especially those design e d to address specific Levy processes, such as VG and 
CGMY, we refer to Hirsa and Madan (2004]), Forsvth et ah ( 2007h . and references therein. 


A.4 Closed-Form Probabilities 

The probability p(t,x ) that an ESO cost surpasses a given threshold x, as defined in (15.31) . can be 
viewed as a European digital option with zero interest rate, computed under P. We summarize the 
corresponding closed-form formulas under the GBM, Merton, and Kou models (see Table [2D • 

(i) Under the GBM model, the ESO cost exceedance probability is given by 

x — x + n(T — t ) 


p(t,x) = e 


- „-A(T-t) 


4>(d), where d = 


a\J T — t 

and 4> is the standard normal c.d.f. 

(ii) When the company stock price follows the Merton jump diffusion, we have 

-A (f-t) e" Q ( T -*)(a(T - t)) y f x -x + p{T - t ) + j 


P (t,x) = c~ A(T -‘T C ' 

3=0 j! \Ja 2 (T -t)+ j a 2 

(iii) In the Kou jump diffusion model, the probability of cost exceedance is given by 


p(t,x) =e 


—Xv 


e (<ry+) 2 v /2 


v n=l k= 1 


-1 


a 

(arj-) 2 v/2 00 


y/m+) k X h-i(x - x - yv; -r)+, -, -arj+y/v) 


av 


+ 


E-nEQnA*V^) k x h-i{x~x- (j,v,-rj- } — ,-ari-y/v) 


a 


tv 

v n =1 


k =1 


av 


, MV — X + X . 

+ 7T0^(--7=- 

ayjv 


where 

n —1 

Pn.fc = ^ 
i=k 
n —1 

Qn,k 

i=k 


i — k 


(A.26) 


(A.27) 


(A.28) 


n A: A AA ( _2±_ ) i-k(_!L ) n-y (1 _ p )"-i, 1 < k < n - 1, P n , n = p n , 
'H+ + V- " 


E 


n — k — l\ /n 
i — k 


V+ 


i)^ , n+ + v-‘ 


v + + v-‘ 

f V- 

y v+ + v- 


y-‘ y _< (l - p)\ 1 < fc < n - 1, = (1 - p)", 
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In{c] d, b, 5) 



E”=o ^ d ) n ~ iHh i(bc -S) + (|f +1 ^e% + &Q(-bc + 5 + 1) 

zr=o fi^Hhiibc -s) + ^_ r+ iv^ e t+^^ bc _ 5 _d b) 


if b > 0, d / 0, 
if b < 0, d < 0, 


with Hh n (x ) = (n!) 1 f£°(t — x) n e t2 / 2 dt, ir n = e av (av) n /n\ and v = T — t. 
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